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Abstract. We present a set of sum rules relating the low-energy light-by-light scattering to 
integrals of 77 fusion cross-sections and use them to study the hadronic contributions. 



1. Introduction 

Some decades ago a general analysis of the forward Compton scattering amplitude allowed 
Baldin [1], Gerasimov [2], Drell and Hearn [3] to establish first sum rules expressing the static 
electromagnetic properties of the nucleon in terms of its total photoabsorption cross sections. 
The Baldin sum rule relates the sum of the electric ue and magnetic 13 m polarizabilities to an 
integral of the unpolarized photoabsorbtion cross section a: 

1 f^dv 

(^e + I3m = ^J -^^{^)^ (1) 

where v is the photon energy in the lab frame. The Gerasimov-Drell-Hearn (GDH) sum rule 
relates the anomalous magnetic moment k of a spin- 1/2 target to an integral of the helicity- 
difference photoabsorption cross section Act = 173/2 — o'i/2 (here subscripts stand for the value 
of total helicity) : 

= - / — Aa(z.), (2) 
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with e the charge and M the mass of the target. 

A few decades later it has been realized [4~6] that the GDH sum rule applies to a photon 
target, in which case the anomalous magnetic moment is zero by Furry 's theorem, and one 
simply has: 

1 f^^dv 



= 
vr 







0-2(1^) - o-o(i^) , (3) 



where a\ is the 77- fusion cross section with the total helicity A, and v is the photon energy in 
collider kinematics. 

More recently, a systematic derivation of sum rules for light-light system has been done [7], 
resulting in sum rules for the low-energy constants of the Euler-Heisenberg Lagrangian. Some 
details of the sum rules for light-by-light (LbL) scattering will be given below, with emphasis 
on perturbative verifications of the sum rules in field theory. The case of charged massive spin-1 
field is especially interesting. 



In the concluding section we shall discuss the way the sum rules can be used to evaluate the 
hadronic contributions to the low-energy interaction of light with light. 



2. Sum rules for light-by-light scattering 

The sum rules in question are based on very general properties of the S-matrix, namely 
Lorentz and crossing symmetries, analyticity, unitarity, as well the gauge symmetry of the 
electromagnetic interaction. To demonstrate this we highlight here the main derivation steps 
for the case of 77 systems. 

Denoting LbL scattering (77 — )■ 77) corresponding Feynman and helicity amplitudes as, 
respectively, M. and M, their relation is: 

where are the photon polarization 4- vectors, A's are the helicities; for real photons traveling 
along the z axis, i.e. q = (0, 0, i/), the polarization vectors are £a(±0 = 2-V2(0, ^A, -i, 0). The 
Mandelstam variables are defined as s = {qi + q-i)^ = Av^, t = {qi — gs)^, u = {qi — 94)^, with qi 
the photon 4-momenta. 

In the forward kinematics, where qs = gi, 54 = q2, and hence t = 0, u = —s, the general 
Lorentz structure of the Feynman amplitude is given by: 

■^/ilM2M3M4 ~ ^{^) 9t^4t^2 9fJ,[it^l ~^ ^{s) 9lJ,4IJ-l9lJ,3H2 

+ C'(s)5M4M3fi'M2Ml ' (5) 

where 5^,^ is the Minkowski metric. Crossing symmetry (under 1 -f-)- 3, or 2 -f-)- 4) means in this 
case 

+ Ciu) 5mim45m3/x2 > (6) 

hence A{—s) = A{s), B{—s) = C{s). As a result there are three independent nonvanishing 
helicity amplitudes: 

M++++(s) = A{s) + C{s), 

M+_+_(s) = A{s) + B{s), (7) 
M++__(s) = B{s) + C{s), 

satisfying the following crossing relations: Af+-+-(s) = M++++(— s), and M+_| (s) = 

M++__(-s)._ 

The principle of (micro-)causality implies that the above functions are analytic functions of 
s everywhere in the complex s plane except along the real axis. For the amplitudes, 

= M++++(s)±M+_+_(s), (8a) 
g{s) = M++__(s), (8b) 

the analyticity infers the following dispersion relations: 
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where j- indicates the principal- value integration. These relations hold as long as the integral 
converges, and otherwise subtractions are needed. Because f^^\—s) = ± f^^\s) and g{—s) = 
g{s), we can express the right-hand side as an integral over positive s only: 
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Re/(-)(.) = -^/d/^^^#^J^. (10b) 



In the physical region (s > 0), the optical theorem relates the imaginary part of these amplitudes 
to the total absorption cross-sections with definite polarization of the initial 77 state: 

Im/(±)(s) = ~[ao{s)±a2is)], (11a) 
lmg{s) = -l[a\\{s)-a^{s)]. (lib) 

Substituting these expressions in the above dispersion relations one obtains: 
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^ -^/d.-.-- ""'^!:-;'-' , (120, 



where a = (ctq + 0"2)/2 = (c7|| + cr±)/2 is the unpolarized total cross section, and Act = (T2 — ctq 
(0 or 2 show the total helicity of the circularly polarized photons, while || or ± show if the linear 
photon polarizations are parallel or perpendicular). 

We next recall that gauge invariance and discrete symmetries constrain the low-energy 
photon-photon interaction to the Euler-Heisenberg form [8], given by the following Lagrangian 
density: 

jO^k = ciiF^uF^^ + C2{F^,F''^f, (13) 

where -F^jy = df^A,y — d^A^, F'^'^ = e^^^^ daAp. Expanding the left-hand side and right-hand 
side of Eq. (12) in powers of s and matching them at each order yields a number of sum rules. 
At 0th order in s we would find 
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= j ds[a\\{s)±a^{s)], (14) 







which cannot work for since the unpolarized cross-section is a positive-definite quantity. 
Empirically a shows a slowly rising behavior at large s and the integral diverges. The assumption 
of an unsubtracted dispersion relation is violated in this case. For the "— " case the sum rule is 
broken too, cf. [5] and references therein. 





Figure 1. Pair production in scalar QED. 



At the first and second orders we find, respectively: 

oo 

= I ds , (15a) 



oo 
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ci ± C2 = — I ds ' ^ — . (15b) 



The first sum rule here is the analog of the GDH sum rule mentioned above, while the sum rules 
for the low-energy constants are unique to the 77 system. Note that according to these sum 
rules the constants ci and C2 are positive definite, in contrast to some previous predictions in 
the literature [9] . 

Most of the above arguments will equally hold for the space-like virtual photons (qf < 0, 
< 0), if written in a variable which reflects under crossing, e.g. v = s — q\ — q2- The EH 
Lagrangian must however be extended by terms containing d^F^^'^ . For the case when at least 
one is real such terms are absent and the above sum rules hold with only a single modification: 
s u = s — q^, where q^ is the other photon virtuality. 

3. Verifications in perturbative QED 

The leading-order cross-sections of 77-fusion in QED are given by trcc-level pair-production 
diagrams (Fig. 1), which can be easily computed and substituted into the right-hand-side of the 
sum rules in Eq. (15). In QED of a scalar and a spinor (Dirac) particle the superconvergence 
sum rule has thus been verified to leading order the fine-structure constant a. For example the 
spin-1/2 pair-production helicity-difference cross section in 7*7 collision is given by 

Aa(T*^->//)(s,g2,0) = , ^"""L 0{s - 4m^) (16) 

(s — q'^y 
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and is plotted in Fig. 2 for three different values of q^. One sees that in all cases the low- and 
high-energy contributions cancel. The fact that 

00 

r ^A.(rW/,( ..0)^^ (17) 
J s-q^ 

is easily verified for any (f' < Am?. 

Substituting the corresponding linearly-polarized cross section into the sum rule for EH low- 
energy constants, we obtain: ci = ^^^^i , C2 = for the scalar case; ci = C2 = 
for the spinor case. This result agrees with the explicit one-loop calculations of low-energy 
light-by- light scattering, see e.g. [10,11]. 
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Figure 2. Helicity-difference cross section of 7*7 — )• // in QED at leading order, for different 
photon virtualities. 



The case of spin-1 QED, describing a charged and massive vector particle, is not as simple. In 
this case interaction with electromagnetic field can be described in terms of three independent 
structure functions. Considering the minimal and linear non- minimal couplings we start with 
the following Lagrangian density: 

£1 = - \f^,^F^'' - ^W*^W>"' - M^W*]^ + ieW*^A''W'' - ieAf,W*W>"' + e^A^W*^W^+ 

(18) 

where and denote the electromagnetic and vector-boson fields respectively; F^j^^ = 
D^Fy — D^F^ and VF^j^ = D^jWy — D^W^ are the corresponding field-strength tensors, with 
= d^ — ieAf^ the covariant derivative, e the charge and M the mass of the vector boson. The 
parameters li and I2 contribute to the magnetic and quadrupole moments as: 

/^ = (l + ^i)^ and g = (/2-/i)^. (19) 

Computing the cross-section for the tree-level 77 — > W~^W^ process in this theory we find 
that the integral on the right-hand-side of the first (super-convergence) sum rule in Eq. (15) 
diverges, unless li = 1 and I2 = 0. Only for the latter choice of the parameters the integral 
converges and is equal to 0, as it should. This choice of parameters is realized in the Standard 
Model, and gives rise to the so-called "natural values" of the electromagnetic moments. 

Computing for the linearly polarized cross sections (for li = 1 and I2 = 0) and using the 
second sum rule we obtain: 

29a2 ^ 27a2 

~ 160M4 ' - Y^^- (20) 

This is in agreement with the one- loop electroweak correction to LbL scattering (Fig. 3), which 
we obtained by the low-energy expansion of the expressions of Bohm and Schuster [12]. It is 
interesting that to obtain this result in electroweak theory one needs to take care of the Higgs 
mechanism as well as the ghosts, while on the side of the sum rule the calculation is much 
simpler: tree-level production of massive vector bosons, and a dispersion integral. 
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Figure 3. One-loop Feynman graphs in the standard electroweak theory. Here W stands for 
the gauge bosons, (f) - Higgs fields and % denotes Fadeev-Popov ghosts. 

4. Hadronic contributions to the LbL scattering 

While the role of the SRs in QED becomes fairly clear in these perturbative calculations, in 
quantum chromodynamics (QCD), in its non-perturbative regime, it is far less obvious. We can 
gain some insight by looking at individual contributions to the 77 hadrons cross sections. 

4.I. Superconvergence sum rule 

The high-energy behavior of Aa is determined by a i-channel exchange of unnatural parity and 
is expected from Regge theory - in the absence of fixed pole singularities - to drop as 1/s or 
faster [13]. Wc therefore expect the sum rule (15a) to converge. The dominant features of the 77 
to multihadron production comes firstly from the Born terms in the 7r'''7r~ (or K^K~) channels, 
which each separately obeys the sum rule. The largest contributions in the hadronic sector are 
thus expected to come from the resonance production: 77 — > M, with M being a meson. It is 
highly nontrivial to see how the sum rule is saturated in this case. 

As two real photons do not couple to a = 1~ or 1+ state due to the Landau- Yang theorem, 
one expects the dominant contribution to come from scalar, pseudoscalar, and tensor mesons. 
One can express the 77 — >■ M cross section for a meson with spin J, mass niMi and total width 
Ttot, using a Breit-Wigner parametrization, in terms of the decay width F^^^ of the meson into 
two photons of total helicity A = 0, 2, as 



T.(A)-p 

rr-^(.) = (2J + l)167r^ * r2 2 

T.(A) 

^ (2J + l)167r2^5(s-m|f), (21) 

niM 

where the last line is obtained in the narrow resonance approximation. For the pseudoscalar 
mesons, which can only contribute to the helicity-zero cross section, the narrow resonance 
approximation is very accurate and allows to quantify their contribution as shown in Tabic 1. 
For the pion, this value is entirely driven by the chiral anomaly, which allows the expression 
of the -K^ contribution to the sum rule as — a2/(47r/^), with /^r = 92.4 MeV the pion decay 
constant. 

To compensate the large negative contribution to the sum rule from pseudoscalar mesons, one 
needs to have an equivalent strength in the helicity-two cross section, (J2- The dominant feature 
of the helicity-two cross section in the resonance region arises from the multiplet of tensor 
mesons /2(1270), 02(1320), and /2(1525). Measurements at various e'^e~ colliders, notably 
recent high statistics measurements by the BELLE Collaboration of the 77 cross sections to 
tt'^'tt" [14], 7r°7r° [15], r^Tr" [16], and K'^K~ [17] channels have allowed accurate confirmation 
of their parameters. As these tensor mesons were also found to be relatively well described by 
Breit-Wigner resonances, we use Eq. (21) to provide a first estimate of their contribution to the 
sum rule. We show the results in Table 2, both in the narrow width approximation and using 
a Breit-Wigner shape, assuming that the tensor mesons pre-dominantly contribute to (T2, as is 
found by the above-mentioned experimental analyses of decay angular distributions. 





m-M 




J ds Aa/s 




[MeV] 


[keV] 


[nb] 




134.98 


(7.8 ± 0.6) X 10-^ 


-195.0 ±15.0 


rj 


547.85 


0.51 ±0.03 


-190.7 ± 11.2 


n' 


957.66 


4.30 ±0.15 


-301.0 ±10.5 


Sum 77, rj' 






-492 ± 22 



Table 1. Sum rule contribution of the lightest pseudoscalar mesons (last column). The 
experimental values of meson masses ium and 27 decay widths r77 are from PDG [18]. 
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Breit-Wigner 




[McV] 


[keV] 


[nb] 


[nb] 


02(1320) 


1318.3 


1.00 ±0.06 


134 ±8 


137 ±8 


/2(1270) 
/^(1525) 


1275.1 


3.03 ±0.35 


448 ± 52 


479 ± 56 


1525 


0.081 ± 0.009 


7±1 


7±1 


Sum /2,/2 






455 ± 53 


486 ± 57 



Table 2. Sum rule contribution of the lowest tensor mesons. We show both results in the 

narrow resonance approximation (4th column) and using a Breit-Wigner parametrization (last 
column). The experimental values of meson masses niM and 27 decay widths r77 are from 
PDG [18]. 

Comparing Tables 1 and 2, we can see that the contribution of the lowest isovector tensor 
meson composed of light quarks, a2(1320), compensates to around 70 % the contribution of the 
vr*^, which is entirely governed by the chiral anomaly. For the isoscalar states composed of light 
quarks, the cancellation is even more remarkable, as the sum of /2(1270) and /2(1525) cancels 
entirely, within the experimental accuracy, the combined contribution of the 77 and rj'. 

Besides the tensor mesons, the subdominant resonance contributions to the 77 total cross 
section arise from the scalar mesons /o/(t(600), /o(980), and ao(980). A rehab le estimate of 
the scalar mesons requires an amplitude analysis of the partial channels, see e.g. [19]. A future 
study will estimate more precisely the scalar meson helicity-zero contribution to the sum rule, 
and elaborate on the cancellation between the tensor mesons and the (pseudo) scalar meson 
contributions in the sum rule of Eq. (15a). Interestingly, when going to the charm sector, the 
sum rule also implies a cancellation between the rjc meson, whose contributions amounts to 
about —15.5 nb, and scalar and tensor cc states. 

In the case of non-zero virtuality of one of the photons, the cross-section of the process 
7*7 — >■ M is defined not by a Breit-Wigner approximation, but rather by a function of Q^, 
so-called transition form-factor (TFF). Application of the superconvergent sum rule to these 
processes can give a useful information of the electromagnetic structure of mesons. Feynman 
amplitude of the process 7*7 — )• M for the case of pseudoscalar meson has the form: 

= ie%,^^q''q'^F^o^.^{Q^), (22) 

where = — g' ^ denots virtuality of the 7* and ^^07*7 is the 7*7 7r° TFF. 

The leading contribution of the single ir^ production to the superconvergent sum rule is found 

as: 

d^^T^ = (^i^.07.7(Q2) j ' . (23) 



The production of a tensor meson contributes both to the hehcity-0 and hehcity-2 amphtudes. 
However the contribution of hehcity-0 cross-section is usuaUy neghgible, while the hehcity-2 
contribution is given by 

r ^'^Tg^ = i^r.m" + ^ (^) • (24) 

Thus, assuming that in the isovector channel the sum rule is saturated by 7r° (A = 0) and 02 
(A = 2) mesons, we obtain at low the following relation between the TFFs: 



F^.r-yiQ') ^ o^T°7*7(Q') ' (25) 



which, as seen from Tables 1 and 2, is at = satisfied to an accuracy of better than 70%. 

4.2. The low- energy- constant sum rules 
The sum rules 



00 

1 r . ^ii(s) 



c, = ^J<^^^, (26) 




C2 



00 



allow one to assess the size of the hadronic contribution to the low-energy LbL scattering. For 
example, the pseudo-scalar meson production in 77* fusion yields: 

(7|| = 0, 

a± = 4Tr^a^\FMjr{Q^)\Us-^Q^)6{s-mlj), (28) 
which lead to ci = and 

C2 = ^ . (29) 



m 



M 



The corresponding numerical results are given in Table 3. This result, however, does not aid 
much in the problem of hadronic LbL contributions to muon anomaly, (g — 2)^, where the main 
effect comes from the LbL scattering at the hadronic scale, see e.g. [20]. 
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C2 [10-^ GeV-*] 







10.8 


V 





0.7 


r{ 





0.4 



Table 3. Contribution of the light pseudoscalar mesons to low-energy LbL scattering. 
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